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Abstract 

Reinforced concrete (RC) shear walls are popular structural systems to resist lateral forces like wind and 
earthquake. The cracking behavior of these critical structural elements is crucial due to its harmful effects on 
structural performance such as serviceability and durability requirements. A technique that shows promise for 
monitoring RC structures is acoustic emission (AE). This paper presents an experimental investigation of fracture 
processes of two large-scale RC shear using AE parameters. In addition, a novel cluster analysis based on the 
k-means is presented to automatically classify the two fundamental cracks modes (i.e., tensile and shear). 
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INTRODUCTION 

Reinforced concrete (RC) shear walls are widely used in conventional building and safety-related nuclear structures. 
They provide much or all of a structure's lateral strength and stiffness to resist earthquake and wind loadings. The 
cracking behavior of these critical structural elements is crucial due to its harmful effects on structural performance 
such as serviceability and durability requirements. In the past two decades, significant efforts have been made 
toward the development of structural health monitoring (SHM) systems in order to reduce life-cycle costs and im- 
prove safety of civil infrastructures. A technique that shows promise for monitoring RC structures is acoustic emis- 
sion (AE) (Aggelis 2011, Carpinteri et al. 2007, Kurz 2006). In general, AE are stress waves caused by sudden strain 
releases due to internal fracture such as concrete cracking (Grosse and M. Ohtsu 2008). Extracting features from 
the AE signal is usually referred to as parameter-based technique. Figure 1 shows a typical AE signal with the most 
common features used during an AE testing (i.e., peak amplitude, rise-time, duration, and count). A threshold is set 
somewhat above the background level and serves as a reference for the waveform features. The AE peak amplitude 
is defined at the maximum amplitude of the signal. The number of times the signal rises and crosses the threshold is 
the count of the AE event. The time period between the rising edge of the first count and the falling edge of the last 
count is the duration of the AE event. Finally, the time period between the rising edge of the first count and the peak 
of the AE event is defined to be the rise time. This paper presents results of an AE monitoring of two large-scale RC 
shear walls (SW1 and SW2) subjected to cyclic loading. Some features were extracted from the AE signals to per- 
form the critical tasks of: 1) crack propagation monitoring and 2) crack classification. In particular, one of the most 
important AE parameters, that is, the b- value has been used to monitor the crack propagation (Farhidzadeh et al. 
2012a, under review) whereas a cluster analysis, based on the &-mean, has been used to perform a crack classifica- 
tion (Farhidzadeh et al. 2012b, under review). This paper is organized as follow: section 2 provides a brief overview 
of the Z?-value analysis and Gaussian filtering. The crack classification procedure based on &-mean analysis is then 
described in section 3 . Section 4 describes the experimental setup followed by the results in section 5 . 
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fc- VALUE ANALYSIS 

The m- value is obtained using the frequency-magnitude distribution data by means of Gutenberg-Richter 
relationship, which is generally used in seismology. The Gutenberg-Richter (1949) formula in terms of AE 
technique is as follow: 



logN(A dB )=a-b(A dB /20) 



(1) 



where N(Acib) is the incremental frequency (i.e. the number of AE events with amplitude greater than the Am), a is 
an empirical constant, b is the Z?-value, and AdB is the peak amplitude of AE event in decibel. The Z?-value analysis 
of AE signals is in general applied to a certain number n of AE signals. Suggested values of n found in literature 
range from 50 to 100 (Shiotani et al. 2001). Several studies have shown that during the process of fracture of an RC 
specimen there is a relationship between the b- value and types of cracks. In general, the b- value increases during 
the formation of micro-cracks (i.e., early stages of damages) and decreases when the macro-cracks begin to localize 
(Colombo et al. 2003). Although Z?-value analysis show promise in monitoring the dynamics of cracks, only few if 
any such methods have been transitioned to field applications. A reason, for this lack of acceptance are the highly 
fluctuations of the b- value caused by the normal operating condition of the structure, which make extremely difficult 
interpreting the results. To overcome this limitation a Gaussian filter is proposed in this work. 



Gaussian Filter 

Gaussian functions have the following properties that make them particularly useful in smoothing filters 
(Lin et. al 1996): (1) The gaussian function is symmetric about the mean, and the weights assigned to signal values 
decrease gradually with distance from the mean. (2) The width of the gaussian function is determined by its spread 
parameter, i.e., the standard deviation. (3) The local extrema (e.g. b- value drop) observed at one standard deviation 
are also observable at the smaller standard deviations and no more local extrema are created as the spread parameter 
increases. Gaussian smoothing is literally the convolution of a gaussian window and a ID vector of data. The 
Gaussian smoothing F(x) of a one-dimentional signal, f(x), is defined as: 



F(x) = f(x)*g(x,a)= jf(fi)g(x-ti 9 <T)dju= jf(ju)-r^exp 



2kg 



2a 2 



dji 



(2) 



where "*" denotes convolution with respect to x, g(x,a) is the gaussian function with the standard deviation a, and ju 
is a dummy variable. These filtering parameters in addition to window span should be properly selected to clarify the 
trend. Example of a local extrema (i.e., b- value drop caused by loading) is shown in Figure 1 . A uniform smoothing 
can eliminate this local drop, depend on window span, but gaussian smoothing will maintain it because it gives more 
weight to window center and less to the tails. 
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Figure 1: Local drop in b- value as a result of loading. 
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Figure 2: Acoustic emission signal with features. 



CRACK CLASSIFICATION 

In the last two decades several efforts have been made to correlate AE features with the cracking mode in RC 
members, which adversely affects structural performance in various ways such as durability and serviceability. Ohno 
and Ohtsu (2010) proposed a technique to perform crack mode classification in concrete structures based on two AE 
parameters, namely: "average frequency (AF)" and the "RA" value, which are defined as (see Figure 2): 



RA = (Rise time) I (Peak Amplitude) 

Average Frequency - (AE ring-down count) I Duration 

Cracks will be categorized as Figure 3 illustrates. However, a defined 
criterion (i.e., the inclined line) on the proportion of the RA value and the 
average frequency for crack classification has not been confirmed yet (Ohno 
and Ohtsu 2010). In this work, a &-means pattern recognition method is 
utilized to classify the AE events in two dominant shear and tensile groups. 

A>mean Clustering 

This method aims at partitioning the data sets into k disjoint subsets 

(clusters) based on minimizing the summation of square distance of each 

data point in a subset to the center of the subset in which it is partitioned 

(Likas et al. 2003). Indeed, several studies based on pattern recognition of 

AE events have determined that it is feasible to distinguish efficiently the 

events associated with various failure mechanisms. Suppose that a data set 

X={xi, X2,... xn}, xn£R d , is available. &-means aims at partitioning this data sets into k disjoint subsets, Ci, . . . , Ck, 

through minimizing the criterion or clustering error as follow 



(3) 



(4) 




RA value 



Figure 3: Signal classification 
based upon JCMS-IIIB5706 code. 



^(w l ,m 2 ,...,m t ) = ^ 1 J M /(^eC / )||x < -m y | 



(5) 



where ra, is the center of/* cluster, Q, and I(P)=1 if P is true and otherwise. This criterion can be minimized 
through iterative procedure summarized in the following steps: 1) Initializing k centers, {mi,. . . , mt}, through 
randomly dividing the data set into k groups and calculating the mean value for each subset. 2) Assigning each 
data point, x n , to a cluster, Q , that its center, mj, has the lowest Euclidean distance to the point among other cluster 
centers. 3) Re-computing the centers for k clusters, {mi,. . . , mt}. 4) If the change in centers in two proceeding steps 
is less than a certain threshold then algorithm is converged and terminate the iteration; otherwise go to step 2. 
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EXPERIMENTAL SETUP 

The test specimens were two large scale RC shear walls, SW1 and SW2, with the thickness of 0.2 m, width of 3.05 
m and height to width ratio of 0.94 and 0.54, respectively. More details about the design of these walls are available 
in (Rocks 2012). The main components of the AE system included an eight-channel high-speed data acquisition 
board (Physical Acoustics Corporation Micro-II PAC) and AEwin software for signal processing and storage. SW1 
was instrumented with eight R15a AE sensors whereas SW2 with four R15a and four R6a. Sensors were attached 
to one face of the wall using hot glue. Preamplifiers were set at 40 dB gain and analog bandpass filters were adjusted 
in the interval of 20 kHz to 400 kHz. Trigger levels of 35 dB and 48 dB were selected to remove the background 
noise for SW1 and SW2, respectively. The sensor layout used for the two specimens is shown in Figure 4. 




Figure 4: RC wall, load actuators, dimensions, and sensor layout for SW1 and SW2. 

Attenuation measurements were carried out by pencil lead breaks made at the center of the specimen, as shown in 
Figure 5. A circular array of AE sensors (R15a from Mistras) were used to estimate the attenuation coefficient in 
different directions, as shown in Figure 5(a). Figure 5(b) shows the topography of the AE amplitudes obtained in the 
SW1 . It can be observed that larger amplitudes were observed in the horizontal direction rather than the vertical one. 
Non-homogeneity and rebar directions of reinforced concrete are the main reason for not having smooth attenuation 
boundaries. The variation range of the attenuation coefficient for the SW1 was between 20 dB/m to 40 dB/m with an 
average of 34.7 dB/m. This average was 29.6 dB/m for the SW2. 




Figure 5: (a)AE sensor layout for attenuation measurements in SW1; 
(b)topography of AE peak amplitude due to pencil lead breaks at the center. 

Table 1 summarizes the coordinates of sensors for loading experiment. The in-plane quasi-static reversed cyclic 
loading protocol consisted of 10 load steps (LS) starting from LSI to LSI for SW1 and 10 load steps starting 
from LS2 to LSI 1 for SW2. Each load step consists of two cycles. The force-displacement hysteresis loops and 
their corresponding backbone (Farhidzadeh et al. 2012a,b; under review) show that onset of nonlinear behavior and 
maximum capacity were observed in LS7 and LS9 respectively, for SW1 whereas in LS8 and LS10 respectively, for 
SW2. 
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Table 1: Sensor coordinates. 








Sensor 


SI 


S2 S3 S4 S5 


S6 


S7 


S8 


Coordinates 
[cm] 


SW1 (23,23) 
SW2 (81,32) 


(157,23) (281,23) (77,103) (157,103) 
(153,32) (225,32) (122,59) (184,59) 


(228.103) 

(81,86) 


(183,156) 
(153,86) 


(50,183) 
(225,86) 



EXPERIMENTAL RESULTS 

To monitor the fracture process of both specimens during the test a Z?-value analysis was carried out as explained 
in section 2. In particular, the Z?-value was calculated using groups of n=70 AE signals. Figure 6 shows the Z?-value 
results obtained considering the entire dataset (i.e. b- values obtained using data collected by all sensors) without the 
filtering process, along with the corresponding probability density function. As can be seen in Figure 6, it is very 
difficult to identify a trend in the unfiltered Z?-value data; this difficulty most likely occurred because of the large 
amount of AE activity (mostly scattering effects) generated by the specimen under the cyclic loadings. A similar 
behavior was observed in SW1. 
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Figure 6: Unfiltered b- values data from SW2. 



Figure 7 shows the results of the gaussian filtering applied to the Z?-value data for the two specimens. A dashed line 
is superimposed to highlight the overall b- value trend. It can be observed that an increasing b- value trend is followed 
by a plateau and then a decreasing trend. 
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Figure 7: Filtered lvalue for (a) SW1; (b) SW2. 

We can infer that micro-cracks formation were dominant prior to LS3 (the increasing trend) and began to localize as 
macro-cracks prior to LS5 (plateau region). From this load step onward, b- values started to decrease (i.e., the macro- 
cracks opening). The decreasing trend began one load step prior to the yielding and reached to a minimum value 
at load step 9. Interestingly, both specimens reached to their maximum strength in the proximity of the load step 9. 
Figure 8 show the maximum and modal (most frequent) crack width in SW1 . A similar crack propagation pattern 
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occurred in SW2. It can be observed that, the maximum crack width was less than 1 mm between LSI and LS6. 
Although numerous cracks were formed in the specimen during these load steps, the crack widths did not change 
significantly. The maximum crack width recorded in the subsequent load steps increased up to 3 mm in LS7 and 
about 4 mm in LS9 as the specimen reached its maximum capacity. It should be noted there is an interest correlation 
between crack widths and b- values, that is, the Z?-value began to decrease at the load step in which the cracks widths 
increased. It can be concluded that a decreasing trend of the b- value may indicate an accumulation of damage in RC 
shear walls. 
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Figure 8: Maximum and modal crack widths, SW1 (Farhidzadeh et. al 2012b, under review). 



Figure 9 and Figure 10 show smoothed b- value analysis results carried out on data recorded by some sensors for 
SW1 and SW2, respectively. The reducing trend and the corresponding load step confirm the conclusion obtained 
considering the data from all sensors. 
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Figure 9: Sensor-based ^-values for SW1 (Farhidzadeh et. al 2012b, under review). 
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Figure 10: Sensor-based ^-values for SW2. 

In concrete structures, tensile cracks appear at low levels of loading while shear cracks appear in severe loading 
conditions. As expected, both diagonal shear and flexural cracks were formed in specimen SW1 due to its 
intermediate reinforcement ratio of 0.67% and relatively large aspect ratio of 0.94. In SW2, a small amount of 
horizontal flexural cracks were generated due to its relatively high reinforcement ratio of 1 .0% and intermediate 
aspect ratio of 0.54. Therefore crack mode classification is more important to be investigated in SW1. 

As discussed earlier, crack mode classification is done by &-mean clustering analysis for SW1 . A moving average 
with window span of 70 hits carried out over the entire RA and AF values in each load step. The related results are 
presented in Figure 11 . It shows that shear cracks become dominant after LS5 and grow quickly afterward. Visual 
inspection about the amount of horizontal and diagonal cracks confirmed such estimation (Farhidzadeh et. al 2012b, 
under review). 
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Figure 11. Crack mode classification for SW1 (a) in LS3; (b) in LS6; 
(c) number of data points associated to tensile or shear mode for the entire experiment 
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CONCLUSIONS 

Reinforced concrete (RC) shear walls are popular gravity and lateral force resisting systems. Failure of a shear 
wall could result in severe damage and even progressive collapse of a concrete structure. This paper described the 
monitoring of the fracture process of two large-scale RC shear walls subjected to cyclic loading using a sparse array 
of acoustic emission (AE) sensors. In particular, a Z?-value analysis was presented to monitor the crack propagation 
during the test. A gaussian filter was proposed to improve the interpretation of the Z?-value data. In addition, a cluster 
analysis based on the &-mean was used to automatically classify the two fundamental crack modes (i.e., tensile and 
shear). 
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